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Abstract 

We conducted radio detection observations at 8.4 GHz for 22 radio-loud broad absorption line (BAL) 
quasars, selected from the Sloan Digital Sky Survey (SDSS) Third Data Release, by a very-long-baseline in- 
terferometry (VLBI) technique. The VLBI instrument we used was developed by the Optically ConnecTed 
Array for VLBI Exploration project (OCTAVE), which is operated as a subarray of the Japanese VLBI 
Network (JVN). We aimed at selecting BAL quasars with nonthermal jets suitable for measuring their 
orientation angles and ages by subsequent detailed VLBI imaging studies to evaluate two controversial 
issues of whether BAL quasars are viewed nearly edge-on, and of whether BAL quasars are in a short-lived 
evolutionary phase of quasar population. We detected 20 out of 22 sources using the OCTAVE base- 
lines, implying brightness temperatures greater than 10 5 K, which presumably come from nonthermal jets. 
Hence, BAL outflows and nonthermal jets can be generated simultaneously in these central engines. We 
also found four inverted-spectrum sources, which are interpreted as Doppler-beamed, pole-on-viewed rela- 
tivistic jet sources or young radio sources: single edge-on geometry cannot describe all BAL quasars. We 
discuss the implications of the OCTAVE observations for investigations for the orientation and evolutionary 
stage of BAL quasars. 

Key words: galaxies: active — galaxies: jets — quasars: absorption lines — radio continuum: galaxies 
— techniques: interferomctric 



1. Introduction 

Broad absorption line (BAL) quasars are a subclass 
of active galactic nuclei (AGNs) with rest-frame ultra- 
violet spectra showing absorption troughs displaced blue- 
ward from the corresponding emission lines in the high- 



ionization transitions of Civ, Siiv, Ny, and Oiv, and 
occasionally in low-ionization transitions, e.g., Mgn and 
Aim (Weymann ct al. 1991). The absorption troughs are 
broader than 2000 km s _1 , sometimes as broad as ~ 0.1c, 
and are presumably due to the intervening components of 
the outflow originating in the activity of central engines. 
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The fact that the most luminous quasars showing BALs 
more frequently (Ganguly et al. 2007) is consistent with 
the strong radiation-prcssure-driven outflows suggested by 
simulation-based studies for accretion phenomena (e.g., 
Proga et al. 2000; Ohsuga 2007). The observed maximum 
velocity of absorption as function of luminosity has an up- 
per envelope, which can be interpreted as the terminal ve- 
locity of radiative driven wind (Ganguly et al. 2007; Laor 
k Brandt 2002). 

The intrinsic percentage of quasars with BALs is ~20% 
(e.g., Hcwett & Foltz 2003). This percentage means 
that the BAL phenomenon takes one of major roles in 
quasars' activities. However, the principal parameter de- 
termining the finding of BAL features is still unknown. In 
the most widely accepted scenario, this percentage rep- 
resents the covering factor of an outflowing BAL wind, 
which is preferentially equatorial, and the BAL features 
can be observed when the accretion disk is almost edge- 
on to the line of sight, based on spectropolarimetric mea- 
surements (Goodrich & Miller 1995; Cohen et al. 1995) 
and a theoretical disk wind model (Murray et al. 1995). 
However, some radio observations provided a counterar- 
gument to this paradigm. Zhou et al. (2006) and Ghosh 
& Punsly (2007) found several BAL quasars with rapid 
radio variability that indicated very high brightness tem- 
peratures, which require Doppler beaming on jets with 
inclinations of less than 35°, i.e., a nearly face-on view of 
the accretion disk. Becker et al. (2000) found that about 
one-third of the radio-detected BAL quasars showed flat 
radio spectra (a > —0.5, S v oc v a ), preferring pole-on jets 
because this geometry tends to make radio sources core- 
dominated by significant Doppler effect only on nuclear 
jets. Also in optical spectropolarimctry, Brothcrton et al. 
(2006) found an electric vector nearly parallel to a large- 
scale jet axis, implying that BAL outflow is not equatorial, 
in a Fanaroff- Riley Class II (FR II) radio galaxy as a BAL 
quasar. Thus, a pole-on outflow would be necessary for 
at least some of the known radio-emitting BAL quasars. 
These results support an alternative proposal that BALs 
are not closely related to inclinations, and may be associ- 
ated with a relatively short-lived (possibly episodic) evo- 
lutionary phase with a large BAL wind-covering fraction 
(e.g., Briggs et al. 1984; Gregg et al. 2000). Gregg et al. 
(2006) pointed out the rarity of FR II/BAL quasars and 
their observed anticorrelation between the balnicity index 
and radio loudness, and suggested that these properties 
can be explained naturally by an evolutionary scheme. 
Montcnegro-Montes et al. (2008) indicated that many 
radio-emitting BAL quasars share several radio proper- 
ties common to young radio sources like Compact Steep 
Spectrum (CSS) or Gigahertz-Peaked Spectrum (GPS) 
sources. Thus, 'inclination angle' and 'evolutionary phase' 
are two of the most important aspects in understanding 
BAL quasars. 

Very-long-basclinc intcrfcromctry (VLBI) instruments 
in radio wavelengths provide exclusive and crucial op- 
portunities to obtain spatial information about AGNs at 
milli-arcsecond (mas) scales by direct measurement, which 
should also be useful for investigating BAL quasars. The 



inclination can be resolved by determining the viewing 
angle of the jet axis, which is supposed to be perpen- 
dicular to the accretion disk. Using the framework of 
Doppler beaming effects, jet axes for many AGNs have 
been estimated by measurements of jet asymmetry (ad- 
vancing speed, brightness, jet length, etc.) by VLBI imag- 
ing. Age as a radio source can be also estimated by 
measurements of its apparent linear size and expanding 
speed, or the age of rclativistic electrons appearing on 
synchrotron spectrum (e.g., Nagai et al. 2006 and refer- 
ences therein). VLBI observations for BAL quasars have 
recently begun to try to study such phenomena. Jiang 
& Wang (2003) observed three BAL quasars with the 
European VLBI Network (EVN) at 1.6 GHz and sug- 
gested that the jet of J1556+3517 was possibly viewed 
from nearly pole-on because of a flat spectrum and un- 
resolved core, while that of J1312+2319 may be far from 
pole-on because of the two-sided structure, and the in- 
clination of J0957+2356 was unclear because of an unre- 
solved steep-spectrum compact source. Liu et al. (2008) 
observed for eight BAL quasar sample, including both 
of LoBALs and HiBALs and both of steep- and flat- 
spectrum sources, with the EVN+Multi-Elcmcnt Radio 
Linked Interferometer Network (MERLIN) at 1.6 GHz. 
High brightness temperatures and linear polarization in 
their core components implied a synchrotron origin for the 
radio emission. No systematic difference was found in the 
radio morphology or polarization properties between their 
limited number of LoBAL/HiBAL or steep/flat-spectrum 
sources. Kunert-Bajraszewska & Marecki (2007) observed 
1045+352 with the US Very Long Baseline Array (VLBA) 
at 1.7, 5, and 8.4 GHz, and found complicated radio mor- 
phology and a projected linear size of only 2.1 kpc, which 
suggests it might be in an early stage, and consistent with 
the evolutionary scenario. These investigations were still 
inconclusive because the detected jet structures, spatial 
resolutions, and frequency coverages were inadequate to 
definitively determine jet properties, and also because of 
still a small number of objects to conclude the natures of 
BAL quasars. 

In this paper, we report our VLBI observations of 
22 BAL quasars at 8.4 GHz by direct measurement at 
the mas resolution, corresponding to the parsec scale at 
the distance to these sources. Our aim is to find BAL 
quasars that have jets suitable for determining their ori- 
entation angles and ages from jet properties for future 
detailed VLBI imaging studies. In Section 2, we describe 
our selection processes, and we present our observations 
and data reduction procedures in Section 3. We present 
the observational results in Section 4, and discuss their 
implications in Section 5. In Section 6, we summarize the 
outcome of our investigation. Throughout this paper, a 
flat cosmology is assumed, with Hq = 70 km s _1 Mpc -1 , 
M = 0.3, and Q A = 0.7 (Spergel et al. 2003). 

2. Sample 

The Sloan Digital Sky Survey (SDSS) Third Data 
Release (DR3; Abazajian et al. 2005) cataloged 



Tabic 1: Optically-selected, radio- flux- limited BAL quasar sample for OCTAVE observation. 
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Col. (1) SDSS source name; Col. (2) Radio counter part from FIRST; Col. (3) Difference between SDSS-FIRST positions; Col. (4) Redshift; Col. (5) Linear scale 
corresponding to 1 mas; Col. (6) BAL subtype, listed in Trump et al. (2006). The code "n" denotes a relatively narrow trough; "HL" denotes a HiBAL in which broad 
(>1000 km s -1 ) low-ionization absorption is also seen; "Hi" denotes a HiBAL-only object, in which broad low-ionization absorption is not seen even though Mgn is within 
the spectral coverage; and "H" denotes a HiBAL in which the Mgn region is not within the spectral coverage of SDSS or has a very low signal-to-noise ratio and so whether 
or not the object is a LoBAL as well as a HiBAL is unknown; Col. (7) i-band absolute magnitude listed in Trump et al. (2006) and calculated using a flat cosmology of 
H = 70 km s" 1 Mpc" 1 , Q M = 0.27, and C A = 0.73 (Spergel et al. 2003); Col. (8) 1.4 GHz peak intensity from FIRST data with a - 5" resolution; Col. (9) 1.4 GHz flux 
density from FIRST data; Col. (10) Radio loudness, the ratio of radio-to-optical flux densities, R* = /5GHz// r (e-g-> Stocke et al. 1992), which were calculated from z, 
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46,420 quasars (Schneider et al. 2005), including 
4784 BAL quasars (Trump et al. 2006). In the catalog of 
the Very Large Array (VLA) 1.4-GHz Faint Images of the 
Radio Sky at Twenty cm (FIRST) survey (Becker et al. 
1995), we searched for radio counterparts with >1 mJy 
located within 10" of these SDSS BAL quasars. This 
search radius was significantly large relative to the posi- 
tion uncertainties of the SDSS and FIRST (~1"), in order 
to rescue even moderately resolved lobe-dominated radio 
sources. Of all SDSS BAL quasars, 91.4% (4374/4784) 
were in the regions of the sky covered by the FIRST sur- 
vey, and we found 492 radio sources corresponding to 
11.2% of SDSS BAL quasars. Most of them were unre- 
solved at a spatial resolution of ~5" of the FIRST, sug- 
gesting intrinsic sizes of less than 1" at the FIRST im- 
age sensitivity (see also Zhou et al. 2006). We selected 
out 23 sources with >100 mJy. Except for three sources, 
these radio sources were found with separations of <1". 
FIRST J004323.8-001548 and FIRST J151005.4+595856 
showed separation angles of 7."6 and 4"3 from their optical 
positions, respectively. The two sources had significantly 
resolved radio structures, which provide the relatively 
large separations; FIRST J004323. 8-001548 had been al- 
ready identified as an FR II-BAL quasar by Gregg et al. 
2006). FIRST J103434.4+592445 was found at 9"9 from 
SDSS J103433.49+592452.3; this was probably a mis- 
identification because an another optical source, SDSS 
J103434. 20+592445. 8 (a candidate of gravitational lcns- 
ing), was a much more likely source for the FIRST source. 
Hence, we removed this source from the list. Finally, 
our optically-selected, radio-flux-limited sample for the 
present VLBI observations consists of 22 sources (Table 1). 

The redshift range is 1.72-3.488 with a median of ~ 2.2 
for the 22 selected BAL quasars. The luminosity and an- 
gular distance are 17456 Mpc and 1705 Mpc at z = 2.2, 
respectively. An angular scale of 1 mas corresponds to a 
linear scale of 8.3 pc at z — 2.2 (see Table 1 for individual 
sources). We determined radio loudness according to the 
definition of the ratio of radio-to-optical flux densities, 
R* = /5GHz// 250Q ^ (e.g., Stocke et al. 1992), which was 
calculated from the redshift z, the FIRST flux density, and 
the SDSS i-band absolute magnitude, assuming a = —0.5 
for the radio and optical spectral indices (see Table 1). 
All of our samples are very radio-loud (logi?* ^ 2.5; 
cf. the extensive study of radio- loudness for FIRST/SDSS 
sources by Ivezic et al. 2002). As this definition may 
somewhat exaggerate the radio loudness of BAL quasars 
with substantial optical reddening. The largest extinction 
At = 0.22 mag is estimated from E(B - V) = 0.114 on 
SDSS J080016. 09+402955.6 in the sample, and exagger- 
ate the radio loudness by log AR* « 0.1. It is also useful to 
consider the rest-frame monochromatic radio luminosity, 
where L 5 qh z > 10 25 W Hz -1 distinguishes radio quasars 
(e.g., Miller et al. 1993). Our samples have L 5 qh z ~ 10 27 - 
10 28 W Hz -1 assuming a = —0.5, are also diagnosed with 
radio-loud objects. 



3. Observations and data analysis 

We observed the selected 22 BAL quasars using 
an optical-fiber-linked real-time VLBI instrument con- 
structed by the Optically ConnecTed Array for VLBI 
Exploration project (OCTAVE; Kawaguchi 2008) and op- 
erated as a subarray of the Japanese VLBI Network (JVN; 
Fujisawa 2008). The OCTAVE connects six radio tele- 
scopes, the Usuda 64 m, Nobcyama 45 m, Kashima 34 m, 
Yamaguchi 32 m, Tsukuba 32 m, and Gifu 11 m, using 
optical fibers at 2.4 Gbps, and offers mJy-level fringe de- 
tection sensitivities. Its high baseline sensitivities are cru- 
cial for efficient exploration of the VLBI-detectability of 
a large number of weak radio populations, such as BAL 
quasars, in preparation for detailed imaging studies. 

Observations were performed in three sessions (Table 2) 
using the Usuda 64 m, Kashima 34 m, Yamaguchi 32 m, 
and Tsukuba 32 m telescopes, which had 8.4-GHz re- 
ceiving systems in the OCTAVE array. Radio signals of 
left-hand circular polarization were received and ampli- 
fied with cooled low noise amplifiers (LNAs) on each an- 
tenna. An analog bandpass filter put through a band- 
width of 512 MHz (8192-8704 MHz), which was sam- 
pled into 2 bits (4 levels) at 1 Gsps using the ADS- 
1000 sampler. The sampled data were transmitted di- 
rectly through the optical fibers to dedicated real-time 
processing correlator units at the National Astronomical 
Observatory of Japan (NAOJ), Mitaka, Tokyo. The 
correlated data were converted into the Flexible Image 
Transport System (FITS) format for analysis using the 
Astronomical Image Processing System (AIPS; Greisen 
2003) software. We were able to obtain all final results 
within an hour after the end of an observation. 

The scan period for each target was 10 minutes includ- 
ing an antenna slew, resulting in on-source integration 
of ^8.0-9.5-minutes. One short scan on a flux calibra- 
tor, which is an almost perfectly point radio source in 
the OCTAVE baselines, was also performed during each 
session (sec Tabic 2). The flux densities of the flux cal- 
ibrators were also measured using a single-dish mode on 
the Yamaguchi 32 m telescope within several days before 
or after the VLBI observations. The flux-calibrator scan 
determined antenna-gain ratios relative to each telescope 
and an absolute amplitude scaling factor at the time of the 
scan. System-noise temperatures were monitored at each 
antenna to correct the time variation of sensitivities. Such 
a method can achieve amplitude calibrations with a dis- 
persion of <5% for antenna-gain ratios and an uncertainty 
of ^10% for an absolute flux scaling, according to past 
Japanese VLBI Network observations at 8.4 GHz (e.g., 
Doi et al. 2006; Doi et al. 2007). The visibility data, which 
had been fringe-fitted and bandpass-calibrated, were av- 
eraged over frequency. 

We have performed non-imaging analysis for these cali- 
brated data. We simply applied a model-fitting to the vis- 
ibility amplitudes of each baseline with a structure model 
of a Gaussian profile. Free parameters of the Gaussian 
were the peak amplitude and full width at half maxi- 
mum (FWHM). The OCTAVE had baselines in the range 
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Date 


Telescope 




Target 




Flux reference 


(1) 
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(3) 
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Col. (1) Observation date (and time); Col. (2) Telescopes used: "U" denotes Usuda 64 m, "Y" denotes Yamaguchi 32 m, 
"T" denotes Tsukuba 32 m, and "K" denotes Kashima 34 m; Col. (3) Truncated name of target source; Col. (4) Flux 
calibrator used for amplitude scaling. 



of 1.5-25 MA, resulting in a sensitivity to sources' FWHM 
sizes ranging from ~2 to 30-100 mas. Because all of the 
baselines are almost in the east-west direction, our one- 
scan snapshot for each target had sensitivity to source's 
brightness profile roughly in one direction only. The cal- 
ibration uncertainty was less than 5% in terms of ampli- 
tude gains relative to each antenna, and baseline sensi- 
tivities of 0.4-1.4 mJy. To avoid amplitude biases due 
to low signal-to-noise ratios, we time-averaged visibilities 
with time intervals sufficient to improve the signal-to-noise 
ratios to more than 10 in the cases of weak sources. 

The results of the model-fitting are listed in Table 3. 
The source sizes of targets detected in more than 
one baseline could be measured straightforwardly. For 
FIRST J004323.8-001548and FIRST J080016.0+402955, 
we could not determine the source sizes because 
no fringe was detected in any baseline. For 
FIRST J021728.6-005226 and FIRST J152821. 6+531030, 
only the shortest baselines could detect fringes, which 
provided lower limits of the source sizes, and also pro- 
vided upper limits with the possibly reliable assumption 
that their total flux densities were equivalent to their 
FIRST peak intensities. For FIRST J100515. 9+480533, 
model-fitting cannot be done because its visibility pro- 
file was quite different from a simple Gaussian. Hence, 
we determined only an upper limit to its source size 
based on the detections with these baselines. For FIRST 
J111914. 3+600457, no model-fitting could be applied be- 
cause data from only one baseline were available, due to 
trouble with the third antenna at that time. We deter- 
mined an upper limit to its source size in the same manner 
as the case of FIRST J100515.9+480533. 

4. Results 

With OCTAVE baselines, we detected 20 of the 
22 FIRST/SDSS BAL quasars. Our OCTAVE observa- 
tions provided baseline detection limits to brightness tem- 
peratures, Tb, of 10 5 -10 6 K for a ~ 9-minute scan period. 
Brightness temperatures were calculated from 

T B = 1.8xl0 9 (l + z) r ^ (1) 



in Kelvin, where z is redshift, S v is the flux density in mJy 
at frequency v in GHz, <f> in mas is the fitted FWHM of the 
source size (cf. Ulvestad et al. 2005). We have derived the 
resolved sizes for 11 sources at mas scales ranging 1.5- 
10 mas, implying Tb = 10 7 1 -10 9 6 K. Five sources were 
unresolved (< 2 mas) with the OCTAVE baselines, im- 
plying Tb ^ 10 8 K. The sensitivity in the determination 
of the brightness temperature for each source depends on 
the projected baseline lengths, signal-to-noise ratio, un- 
certainty in amplitude calibration, and any discrepancy 
between the two-dimensional intrinsic source profile and 
the one-dimensional fitting model (see Section 3). For 
FIRST J100515.9+480533, the visibility profile cannot be 
represented by a simple Gaussian, which may suggest the 
existence of multiple components with comparable fluxes 
at mas scale. A complicated jet structure was inferred in 
this BAL quasar. 

We also show two-point spectral indices between the 
FIRST 1.4-GHz peak intensities and the maximum corre- 
lated flux densities of the OCTAVE at 8.4 GHz (Table 3). 
We used the peak intensities rather than total flux densi- 
ties at 1.4 GHz to extract emitting components as close 
to the nuclei as possible. The FIRST data were obtained 
with a ^"-resolution, which is much larger than that of 
the OCTAVE; thus, the resolution effect could occur in de- 
riving the spectral indices. Note that the derived spectral 
indices provide only lower limits, if any of the radio emit- 
ting components in a source lies outside a region within 
~10 mas. For this reason, we can be confident of the de- 
tections of inverted spectra (a > 0) . On the other hand, we 
cannot determine whether steep spectra (a < 0) were the 
results from an intrinsic nature or a resolution effect. In 
the 22 targets, we found inverted spectra in four sources. 
Note that the derivations of spectral indices could also be 
affected by flux variability, as ~ 10 years elapsed between 
the FIRST and OCTAVE observations. We discuss this 
possibility later. 

5. Discussion 

5.1. Coexistence of nonthermal jet and BAL outflow 

We detected 20 radio sources using OCTAVE base- 
lines, which assured brightness temperatures of greater 



G 



Doi et al. 

Table 3: Results of OCTAVE observations for 22 BAL quasars. 
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Table 3: (continued) 
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Col. (1) Truncated name of target source; Col. (2) Correlated flux density; Col. (3) Projected baseline length; Col. (4) Fitted flux density 
of a Gaussian profile (see Sections 3 and 4); Col. (5) Fitted FWHM size of a Gaussian profile; Col. (6) Fitted FWHM in pc; Col. (7) 
Brightness temperature calculated using Eq. (1); Col. (8) FIRST 1.4 GHz peak intensity, the same as in Table 1; Col. (9) Spectral index 
between 1.4 and 8.4 GHz, calculated from the FIRST peak intensity and the maximum correlated flux density. 
* Undetected at all baselines, indicating a la upper limit of the shortest baseline, 
t Detected only at the shortest baseline. 

t Complex visibility profile; cannot be fitted with a Gaussian. 
§ Only a single baseline observation. 

than 10 5 -10 6 K. Although the brightness temperature 
of a radio supernova at a very early stage could exceed 
T B = 10 6 K (e.g., Bietcnholz et al. 2001), its expected flux 
density would be far from sufficient for fringe detection 
for such distant quasars in our sample (cf. van Dyk et 
al. 1993). The brightness temperatures of stellar compo- 
nents of luminous starbursts are < 10 5 K (cf. Lonsdale 
et al. 1993), which are less than those of the detected ra- 
dio counterparts of BAL quasars. Also, in terms of radio 
luminosity, even the most radio-luminous starbursts show 
up to only ~ 10 24 W Hz" 1 (e.g., Smith et al. 1998), in con- 
trast to ~ 10 27 -10 28 W Hz -1 for our BAL quasar sample. 
Therefore, the OCTAVE-detected radio emissions cannot 
be accounted for by any stellar origin. We conclude that 
the OCTAVE-detected radio emissions of the BAL quasars 
originate in nonthermal jets from AGN activity, as in the 
cases of other AGN radio sources. VLBI observations 
that previously conducted (Jiang & Wang 2003; Kunert- 
Bajraszewska & Marccki 2007; Liu et al. 2008) also sup- 
port the existence of nonthermal jets in BAL quasars. 

The fairly high VLBI-detection rate (20/22) is evidence 
that BAL outflows, which are inferred from broad troughs 
in UV spectra, can coexist with nonthermal jets in radio- 



loud BAL quasars. This indicates that the accretion disks 
of BAL quasars can generate radiation-pressure driven 
strong outflows and magnetic-driven strong jets simulta- 
neously. It is important to investigate the relationship 
between the properties of BAL features and radio emis- 
sions (cf. Ghosh & Punsly 2007) to understand accretion 
phenomena in quasars. 

5.2. Jet properties of observed BAL quasars — inverted- 
spectrum sources 

We found four inverted-spectrum (a > 0) sources 
(Column (9) in Table 3). Since optically-thin nonther- 
mal synchrotron emission should show steep spectrum 
of a ^ —0.5, the observed inverted spectra should result 
from a lower-frequency absorption mechanism if a non- 
simultaneous spectral index is not affected by flux vari- 
ability. Absorbed components should dominate the total 
spectrum for an inverted spectrum to be observed even us- 
ing the FIRST beam width. Hence, an inverted spectrum 
suggests three possibilities: (1) Doppler beaming effect on 
jets, (2) a young (compact) radio source, or (3) artificially 
made due to flux variability, as follows. 

Doppler boosting can apparently enhance only 
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optically-thick nuclear components beyond extended jets 
that would have been decelerated and optically-thin: 
the nuclear components tend to be synchrotron self- 
absorbed because of high-brightness temperatures, and 
show inverted spectrum at lower frequencies. An ade- 
quate Dopplcr beaming effect requires jets that are nearly 
aligned with our line-of-sight; this implies a face-on viewed 
accretion disk, which is inconsistent with the widely ac- 
cepted scheme of equatorial BAL outflows. The presence 
of Doppler beaming effect has already been inferred from 
rapid flux variation between the NRAO VLA Sky Survey 
(NVSS; Condon ct al. 1998) and the FIRST in several 
BAL quasars (Zhou et al. 2006; Ghosh & Punsly 2007). 
It is important to confirm the Doppler beaming in also by 
VLBI imaging in a future. 

Young radio sources also can make their spectra 
nuclear-dominated, because extended jets would not yet 
been developed. Radio sources of ~ 100-1000 pc or 
less could show peaked spectra, resulting inverted spec- 
tra possibly at <~ 1 GHz; the spectral-peak frequency is 
roughly determined by the linear size of radio structure 
(e.g., Snellen et al. 2000). The explanations of the rela- 
tion between linear size and spectral-peak frequency have 
been suggested using synchrotron self-absorption (O'Dea 
& Baum 1997) and free-free absorption (Bickncll et al. 
1997). Montcncgro-Montes et al. (2008) presented many 
radio sources of BAL quasars showing such peaked spec- 
tra. The ages could be determined from the linear size 
and expanding speed of radio structures; it is important 
to measure the expanding speed by VLBI-imaging mon- 
itor in a future. For example, in a typical VLBI scale, 
an apparent expanding rate of 0.1 mas yr -1 and a linear 
size of two-sided structure of 100 pc leads to an estimated 
age of 500 yr. This situation can be adapted to the evolu- 
tionary scenario of BAL quasars in a relatively short-lived 
phase (e.g., Briggs et al. 1984; Gregg et al. 2000, and see 
also Gregg et al. 2006). 

Inverted spectra artificially-made due to flux variability 
between the observations of the OCTAVE at 8.4 GHz and 
the VLA at 1.4 GHz (~ 10 yr) cannot be ruled out. A 
fraction of BAL quasars are significantly variable in the 
radio band (Zhou et al. 2006; Ghosh & Punsly 2007). A 
flux variability of >10% in 10 yr at 8.4 GHz for a 100- 
mJy source at z = 2 implies Tb > 10 113 K, which would 
be above the inverse-Compton limit (Readhead 1994) and 
require Doppler beaming effect. A flux variability of 10% 
corresponds to a change of spectral index of only ~ 0.05. 
That means that even if the observed spectra had been ar- 
tificially made, we have the same conclusion that the four 
inverted-spectrum sources are possibly Doppler-boosted. 
As a result, we conclude that the inverted-spectrum sources 
are interpreted as Doppler -beamed, pole- on-viewed rela- 
tivistic jet sources or young radio sources such as CSSs 
and GPSs: single edge- on geometry cannot descibe all 
BAL quasars. 

We also check flux variation between NVSS and FIRST 
at 1.4 GHz for our BAL quasar sample. Using the same 
manner as Ghosh & Punsly (2007), we found two sources, 
FIRST J075628. 2+371455 and FIRST J122836.9-030438, 



with significant flux variation of 4.8er and 3.2er; brightness 
temperatures were derived to be 10 15 3 K and 10 13 2 K, re- 
spectively. Both of the radio sources were very compact 
(2 mas or less) in OCTAVE observations, and the latter 
showed an inverted spectral index (a = +2.0). These situ- 
ations are consistent with the presence of highly Doppler 
boosting on pole-on jets. It is important to confirm them 
also by VLBI imaging in a future. 

5.3. Jet properties of observed BAL quasars — steep- 
spectrum sources 

The majority of the detected BAL quasars showed steep 
(a < 0) radio spectra. It was unclear whether the weaker 
radio flux densities were observed at 8.4 GHz because ex- 
tended structures were resolved out, or were due to intrin- 
sically steep spectra on compact components. Hence, we 
do not discuss further in detail for the OCTAVE results. 
We have already stressed the importance of VLBI imag- 
ing observations for inverted-spectrum sources, and steep- 
spectrum sources are also needed to be VLBI-observed 
because they are the majority of BAL quasars and are 
also compact in most cases (Becker et al. 2000). The 
correlated flux densities in the OCTAVE baselines were 
not much smaller than the 10 _1 -times the VLA peak in- 
tensities in most sources, despite the fact that the beam 
area of the OCTAVE was - 10~ 5 -timcs that of VLA. This 
indicates that radio-emitting origins considerably concen- 
trated in parsec-scale components, which should be inves- 
tigated using VLBI. As the evolutionary scenario, BAL 
quasars might associate young radio sources (Briggs et 
al. 1984; Gregg et al. 2000; Gregg et al. 2006), which 
should be optically-thin small radio lobes seen in com- 
pact steep spectrum (CSS; O'Dea 1998 for a review) ob- 
jects. Many CSS sources have been revealed as young 
(< 10 5 yr) radio galaxies by VLBI-imaging studies (e.g., 
Murgia 2003; Nagai et al. 2006). VLBI-imaging studies 
may provide evidence supporting the evolutionary sce- 
nario for BAL quasars, if the steep-spectrum radio sources 
of BAL quasars are CSS objects (Kunert-Bajraszewska & 
Marccki 2007). 

5.4- Implications of OCTAVE observations 

Our OCTAVE observations have many implications for 
the study of BAL quasars. VLBI images of only four 
BAL quasars have been published before the OCTAVE 
observations (Jiang & Wang 2003; Kunert-Bajraszewska 
& Marecki 2007). Our OCTAVE observations have dra- 
matically increased the number of VLBI-detected BAL 
quasars, and have established that this AGN subclass 
includes a non-trivial number of radio-loud objects that 
can be directly imaged at mas resolutions (correspond- 
ing to parsec scales) as well as objects in the other 
AGN classes. The orientation angle and the age as ra- 
dio sources of BAL quasars should be determined from 
jet properties in subsequent detailed VLBI imaging stud- 
ies. The OCTAVE strongly recommends that multi- 
epoch and multi-frequency (and polarimetric) VLBI ob- 
servations should be performed for these sources. Multi- 
epoch imaging can measure the proper motion of jets or 



VLBI Detections of Radio-Loud BAL Quasars 



9 



lobes to estimate the inclination or kinematic age. Multi- 
frequency imaging can measure the spectral-index profiles 
along approaching and receding jets or lobes, which can 
discriminate between Doppler-boostcd self-absorbed jets 
and free-free absorbed jets obscured by thermal plasma. 
BAL outflows could be free-free absorber along nonther- 
mal jets; the projected one-dimensional profile of thermal 
BAL outflows could be studied in mas resolutions using 
VLBIs at multi-frequency at by measuring opacity gra- 
dient along the jets. The free-free absorber also could 
occur Faraday rotation to the vector of polarization axis 
toward jets, which is an another powerful tool to investi- 
gate the spatial profile of thermal BAL outflows. In multi- 
frequency VLBI observations, such processes should offer 
exclusive probes to the parsec-scale profile of BAL out- 
flows. On the basis of the results of the OCTAVE obser- 
vations, we arc observing some of the OCTAVE sample 
by multi-frequency VLBI imaging. It will be important 
to compare the results of inclination measurements based 
on the VLBIs and optical spectropolarimetry, and to in- 
vestigate the relation among UV/optical spectra, pc-scale 
geometry, and ages. 

6. Summary 

Two of the most important questions in understanding 
BAL quasars are: (1) Arc they viewed at nearly edge- 
on? (2) Are they in a short-lived evolutionary phase? 
Our OCTAVE observations have increased the number of 
VLBI-dctected BAL quasars, which offers more chance 
to determine their orientations and ages by subsequent 
VLBI-imaging studies to conclude the two pictures. We 
detected 20 out of the radio-brightest 22 sources selected 
from the counterparts of SDSS BAL quasars. We con- 
cluded that nonthermal pc-scale jets and thermal BAL 
outflows can coexist in these radio-loud BAL quasars 
simultaneously. BAL quasars have become targets of 
VLBIs to be revealed in pc scale by direct imaging, 
as in the cases of other AGN classes. We also found 
four inverted-spectrum sources, which are interpreted as 
Dopplcr-beamed, polc-on- viewed relativistic jet sources or 
young radio sources: single edge-on geometry cannot de- 
scribe all BAL quasars. 

We thank M. S. Brotherton for very useful comments 
that improved the presentation of the paper. We also 
thank K. Asada for very useful comments. This work 
was partially supported by a Grant-in-Aid for Young 
Scientists (B; 18740107, A. D.), a Grant-in-Aid for 
Scientific Research (C; 21540250, A. D.) and a Grant-in- 
Aid for Scientific Research (C; 20540233, K. W.) from the 
Japanese Ministry of Education, Culture, Sports, Science, 
and Technology (MEXT). We are grateful to all the staffs 
and students involved in the development and opera- 
tion of the OCTAVE. The OCTAVE project has been 
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of Radio Astrometry (VERA) project in the National 
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the National Institutes of Natural Sciences (NINS). The 



optical-fiber networks for OCTAVE have been provided 
by the GALAXY project 1 supported by NTT Corporation 
(Uose ct al. 2002), the Japan Gigabit Network-2 (JGN2) 
project operated by the National Institute of Information 
and Communications Technology (NICT), and the Science 
Information NETwork 3 (SINET3) operated by the 
National Institute of Informatics (Nil). The OCTAVE ar- 
ray consists of six contributed antennas: the Usuda 64 m 
of the Japan Aerospace Exploration Agency (JAXA), 
the Kashima 34 m of the NICT, the Nobeyama 45 m 
of the Nobeyama Radio Observatory (NRO) in the 
NAOJ, the Tsukuba 32 m of the Geographical Survey 
Institute (GSI), the Yamaguchi 32 m (operated by 
Yamaguchi University) of the NAOJ, and the Gifu 
11 m of Gifu University. Time allocation and array 
operation for OCTAVE are being carried out under 
the framework of the Japanese VLBI Network (JVN) 
project, which is led by the NAOJ, Hokkaido University, 
Gifu University, Yamaguchi University, and Kagoshima 
University, in cooperation with the Institute of Space 
and Astronautical Science (ISAS)/JAXA, GSI and 
NICT. We used the US National Aeronautics and 
Space Administration's (NASA's) Astrophysics Data 
System (ADS) abstract service, the NASA/IPAC 
Extragalactic Database (NED), which is operated by the 
Jet Propulsion Laboratory (JPL). In addition, we used 
the Astronomical Image Processing System (AIPS) soft- 
ware developed at the US National Radio Astronomy 
Observatory (NRAO), a facility of the National Science 
Foundation operated under cooperative agreement by 
Associated Universities, Inc. 
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